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Abstract 

Motivated by recent experimental results, we use a factorization approach to study the three- 
body B D^*^ K'^ decay modes. Two mechanisms are proposed for kaon pair production: 
current-produced (from vacuum) and transition (from B meson). The B^ D^*^^ decay 
is governed solely by the current-produced mechanism. As the kaon pair can be produced only 
by the vector current, the matrix element can be extracted from e~^e~ — > KK processes via 
isospin relations. The decay rates obtained this way are in good agreement with experiment. Both 
current-produced and transition processes contribute to B~ D^*^'^ K~ decays. By using QCD 
counting rules and the measured B^ D^*^^ K~ decay rates, the measured decay spectra can 
be understood. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



The Belle Collaboration reported recently the first observation oi B ^ [)(*) j{- 
decays [0, with branching fractions at the level of 10~^ — 10^^. Angular analysis of the 
K~K^*^^ subsystem reveals that K~K^ and K~K*^ meson pairs are dominantly = 1~ 
and = 1^, respectively. While there is no sign of decay via resonance for the K~K^ pair, 
data suggest a dominant ai(1260) resonance contribution in the production of the K^K*^ 
pair. Both the K'K*^ and the K'K^ pair mass spectra of the D^*^K-K*° and D^K~K^ 
modes show a maximum near threshold. 

These processes can be described by conventional b cud diagrams with additional ss 
pair creation. For example, in the £){*)+ K^*^^ case, the spectator quark (J) ends 

up in the meson, while the ud pair ends up in the kaon pair. It is similar to the 

observed three-body baryonic mode B^ D*^pn 0], where the nucleon-antinucleon pair 
replaces the K~K^^*^ pair in the above. A generalized factorization approach has been 
applied to study this three-body baryonic mode, where the amplitude is factorized into a 
current-produced pn part and a B^ ^ D*^ transition part. 

A crucial ingredient in the factorization approach to the B^ D*^pn mode is the 
knowledge of the time-like nucleon form factors. Although the time-like axial nucleon form 
factor data is still not available hence making the axial current contribution incalculable, the 
vector current induced nucleon form factors can be obtained via isospin rotation from their 
EM counterparts, where data is quite abundant. By utilizing nucleon EM form factor data, 
it was shown that the vector current contribution can account for up to 60% of the observed 
B^ — s> D*^pn rate. The predicted pn mass spectrum shows threshold enhancement effect as 
expected [Q . This can be seen as rooted in the near-threshold behavior of the nucleon form 
factors, whose appearance can be traced back to the application of factorization and QCD 
counting rule, which has been confirmed in the nucleon EM data. The total rate might be 
fully understood once the axial nucleon form factor becomes available. Alternatively, it was 
proposed that one can extract the axial nucleon form factor from future B'^ D*^pn data. 
A factorization and pole model approach has recently been employed to study B ^ D pn, 
D pp modes ||^, and some information on axial form factor is extracted. 

With success in the three-body baryonic modes, and the encouragingly similar threshold 
enhancement behavior in the B [)(*) K~ K^*^'^ modes [|I|, we apply the factorization ap- 
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proach to study these three-body decays. For the B 



D'^*^+K-K° modes, the fact that 



the K^K^ pair is observed only in the = 1~ state aheady supports the factorization 
picture, since only the vector current can produce the kaon pair under factorization. With 
no axial current contribution, the amplitude can be predicted by using kaon EM form factors 
through isospin relations. This is in contrast with the B D*~pn case where the axial 
current also contributes. Thus, with no tuning of parameters, the B^ D^*^~^ K~ modes 
provide a useful testing ground of the factorization approach. 

On the other hand, for the B^ D^*^^ modes, the kaon pair can also be produced 
from a B meson by a current induced transition. The situation is similar to B ^ pp 
transitions in the B ppK case 0. Since the B to kaon pair transition form factors are 
not known, we use a parametrization motivated by QCD counting rules, and determine these 
parameters from total decay rates. In other words, our approach is less predictive for these 
modes compared to B^ D^*^^ K~ . However, the predicted decay spectra are closely 
related to the QCD counting rules, which can be tested experimentally. 

In this work we shall concentrate on the B D^*'>K~K^ modes, while making only some 
comments on the B D^*^ K~ K*^ modes. The experimental data indicate that the K~K*^ 
subsystem is in a 1+ state and dominated by ai(1260) resonance, hence originate from the 
axial current. But, unlike the vector current case for K^K^ production, we do not have 
independent information on the axial form factors, so we have no control over these modes. 

This paper is organized as follows: in the next section we lay down the formalism of 
the factorization approach. We show how to extract kaon form factors from EM data in 
Sec. Ill, and parametrize the transition form factor. In Sec. |V| we show the results of our 
calculation, and in the last section we make some discussions before conclusion is drawn. 

II. FACTORIZATION FORMALISM 

The relevant effective Hamiltonian for the 6 — > c transition is 



where Ci(/i) are the Wilson coefficients and Vcb and Vud are the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements. The four-quark operators Oi are products of two V — A 
currents, i.e. Ol = {cb)v-A {du)v-A and O2 = {cu)v-a- 




v,,v:,[c^{^l)Ol{^^) + cMO^,{^l)], 



(1) 
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With the factorization ansatz, the decay amphtudes for B D^*^K are given by 

A{D'^*^'K-K') = ^K6V;*,[ai(Z}W+|(c6)^_^|5°)(ir-ir°|(Jn)v^_A|0) 

+a2{K-K'\{db)v-A\B-){D'^*'>'\{cu)v-A\0)], (2) 

where the effective coefficients are expressed as ai = ci + C2/3 and a2 = C2 + ci/3 if naive 
factorization is used. The factorized amphtudes consist of products of two matrix elements: 
the case of B to D transition times current produced kaon pair is caUed "current-produced" 
(denoted as J'), while the case of B to kaon pair transition times current produced D is 
called "transition" (denoted as T). The two cases are depicted in Fig. Note that the 
BO D^*^+K~K° decay can only be current-produced, while B D^*^^K receive 
both contributions. 

The {D^*^\V-A\B) and {D^*^V - A\0) matrix elements are familiar and can be param- 
eterized in the standard way. We shall adopt the Bauer-Stech-Wirbel (BSW) and the 
Melikhov-Stech (MS) models for comparison. The matrix elements {KK\V — A\0) and 
{KK\V — A\B) are less familiar. They are parametrized as ^ 

{KMK2{P2) \V>^\ 0) = {pi -p2rF^^{q'), (3) 

{K,{p,)K2{p,)\{V - A)^\B{pb)) = iw4q')iP2~Pi),+ hiq')e,,^pp'i,q'^{p,-p,)^, (4) 

where q = Pi + P2, and we have dropped ml — m\ dependent terms in Eq. (^ by assuming 
isospin symmetry. Since {K~ K^\{du)A\^) = from Lorentz covariance and parity, only 
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the vector current contributes to Eq. (0), which explains the experimental observation of 
= 1~ for K~K'^ pair in D^*^^ K~ . Data also show that the kaon pair for the 

B- D^*^^K-K° modes is also in a 1 configuration Since the transition amplitudes 
now also contributes, we have used this experimental fact to drop the {ps — q)^ and 
dependent terms in Eq. (^), as they would lead to other quantum numbers for the kaon pair. 
This greatly simplifies the work. 

We now show how to use data on time-like kaon EM form factors and an isospin relation 
to obtain F^^{q^). We then give a simple parameterization on 5 — > KK transition form 
factors motivated by QCD counting rules. 



III. CURRENT-PRODUCED AND TRANSITION KK FORM FACTORS 



A. Isospin Relation and Kaon Electromagnetic Form Factor 

The D^*^^ modes contain only current-produced amplitudes J^, as can be seen from 
Eq. (^) and depicted in Fig. Since only vector current contributes, we can use kaon EM 
form factor data to obtain the weak vector form factor via the isospin relation 



(5) 



where , F^o are the EM form factors of the charged and neutral kaons, respectively. 

The kaon EM form factors have been measured for both space-like and time-like re- 
gions [p!0| , pj] , |13[, where processes e^e^ K^K~, KlKs provide the time-like data. 
The time-like \Fx+\ and |-Fft:o| form factor data are given in Fig. ^ in the energy region 
Mk+k-^ MklKs = 1 ~ 3 GeV. The structure is complicated in the 1 ~ 2.1 GeV range, re- 
vealing both resonant as well as non- resonant contributions. A sharp 0(1020) peak is shown 
in the insets. The form factors drop quickly above 1.02 GeV, but a slower damping takes 
over for larger Mkk, and one must include p, u, (p and their higher resonances in modeling 
the form factors. 

The asymptotic behavior in Mkk is characteristic of power-law fall off, and seems to 



satisfy a stringent asymptotic constraint Jill from perturbative QCD (PQCD), 

-1 



Fxit) 



(6) 
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FIG. 2: Time-like (left) and |i^/^o| (right) form factor data, where the inset is for (j) region. 

They are fitted by Eqs. (|^), (^), respectively. 

where t = and A ~ 0.3 GeV is the QCD scale parameter. The 1/t power reflects the need 
for a hard gluon to redistribute large momentum transfer. 

Since our aim is just to parametrize and fit the form factor data, we express the EM form 
factors by a phenomenological model that combines the resonant and PQCD terms 

-1 



FK+{t) 
FK^t) 



t — m?- + inijTjit) 



Xi X2 



In 



t 

A2 
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t — m'j + inijTjit) 



In 



A2 
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(7) 
(8) 



where one sums over appropriate meson poles and Ij is the isospin of the jth meson. We 
put very few asymptotic PQCD terms because the large t data is sparse. The Coulomb 
factor C{t) [|15] accounts for soft photon exchange that can enhance the cross section by a 
few percent for low and is needed for the fit. But in obtaining F^^ from Fk+ via Eq. (^, 
we divide out this effect since there is no exchange of soft photons for K~K^. 



B. Fitted Kaon Weak Vector Form Factor 

A major difference with our earlier studies of baryonic modes 0] is the presence of the 
resonance part, which turns out to be important. To be able to account for the rich structure 
in the experimental data in mK-KO,iTLKi^,Ks ^ 2 GeV, we take the eight vector mesons 
p(770), Lj(782), 0(1020), cj(1420), p(1450), u;(1650), 0(1680), p(1700), with their masses rrij 
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kept at the experimental values [T^. The widths Tj{t) for j = p, u, (p are given by |TT 



snir, 



2s 



t - Ami 
— 4m^ 

t - Am\+ 



3/2 



r^(t) = r. 



3/2 



3/2 



ml - Am\o 



(9) 



where Tp^^^ are the full widths |16]. For the other higher-mass vector mesons, we simply 



take their experimental width values 



We should stress that our phenomenological model is devised to account for EM form 
factor data in the region of our interest. No attempt is made for deeper theoretical un- 
derstanding, nor for data outside the t-region of interest. We therefore do not need all the 
0-region data even though more precise measurements have been obtained, as one would 
need more sophisticated treatment of T^{t) (including a — Svr phase space term) |T^ 
which would complicate our parameterization. Actually, our neglect of such data causes 
no harm since (j) is an isoscalar. While it is needed to account for Fk+ and -FVo, its effect 
cancels in the weak form factor Ff"^{q^)^ as required by isospin symmetry. 

The coefficients Cj are treated as free real parameters though by some physical conditions 
they are not all independent. In the VMD framework, Cp, c^, are expressed as Cp = 
dp-yOpRK, = Qui-yQuiKK and C0 = g^g^iKK- Usiug experimental values of the electronic 



widths [|T6| for the coupling constants gy-y and assuming SU(3) relations with ideal mixing, 
namely gpKK = g<t>KK/ = gujKK, we have ITT 



1:1:1, 
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(10) 



For other higher-mass vector mesons, Cj are free from the above constraint. 

We take Eqs. (|^ and (^ to make a phenomenological fit to the experimental data of the 
charged [0, ^ and neutral kaon form factors [1^. The best fit values are obtained by 
finding the minimum of the function, = X+ + Xo' where 

\FK+4t,)\f 



(11) 



where ZiS are the measured absolute values of the time-like kaon EM form factors and (TjS 
are the error-bars. Note that we fit the absolute values of the kaon form factors since these 
are what can be obtained by experiment. Also, due to the common resonance part in both 
Fx+ and Fj^o, we search for the minimum of as a combination of these two form factors. 
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FIG. 3: The kaon weak vector form factor F^^{t). 



We find the best fit values (in unit of GeV^) 



Cp = 0.363, 



Cp(i450) = 7.98 X 10-3, 



Ca;(1420) = -7.64 X 10 ^ C<^(1650) = "0.116, 



Cp(1700) = 1.71 X 10 ^ 
c<^(i680) = -2.0 X 10"^ 



(12) 



and 



xi = 3.26 GeV^ 



X2 



-5.02 GeV^; yi = -0.47 GeV^ 



(13) 



Note that the ratios of c^, with Cp are already fixed by Eq. ([T^), and the best fit value for 
Cp is consistent with Ref. [|rT|. As can be seen from Fig. 0, the fit is reasonable (x^/n.d.f = 
194/130 ~ 1.5) for both low and high energies. 

Using Eq. (^, we give the kaon weak vector form factor Ff^^ in Fig. ^, with the overall 



factor ^/C{t) in Fk+ removed. Contributions from poles of J = mesons cancel out, and 
one is left with only the p(770), p(1450), p(1700) and PQCD contributions. Since Cp(i45o), 
Cp(i7oo) are small, 

pKK 

is modeled by the p(770) and PQCD parts hence quite smooth. 
Under factorization assumption, this could be the reason behind the absence of structures 
in the K~K^ pair spectrum for the 5° D^*''^K'K^ decays 



C. Ansatz for B — > K Transition Form Factors 



In addition to the current-produced matrix element, we also need the transition matrix 
element {K' K^\{dh)v-A\B~) for B- -> D^*^^K-K^ decay amplitudes. Two out of four 
possible transition form factors are eliminated by the experimental observation of = 1~ 
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for the kaon pair |T|. There is no experimental data on the remaining form factors, denoted 
w^{q'^) and /i(g^) in Eq. (g). 

It is interesting to study the asymptotic behavior of these form factors from QCD counting 
rules. To produce a kaon pair with large invariant mass from a decaying B meson, at least 
two hard gluon exchanges are needed: one creating the ss pair in K^K^, the other kicking 
the spectator to catch up with the energetic s quark to form the K meson. This gives rise 
to a asymptotic behavior, where t is the K~K^ invariant mass squared. 

Resonant contributions such as from an intermediate p pole should in principle be con- 
sidered for the K^K^ transition. This seems to be supported by the dominant p 
contribution in F^^{q^) as we have just discussed. However, we would need further poles 
such as p(1450), p(1700) to account for the asymptotic behavior implied by QCD count- 
ing rules. Our experience with F^^ does not help since these resonances are unimportant 
there, while we lack other independent experimental information. But there is as yet no 
clear sign of resonances in the kaon pair spectrum. Because of this, we shall use a very 
simple parameterization solely motivated from QCD counting rules, i.e. 

^-W = ^, Kt) = ^, (14) 

where c^^h are free parameters to be fitted by data. 
IV. RESULTS 

We use the central values of the effective coefficients = 0.91 ± 0.08 ± 0.07 (0.86 ± 

0.21 ± 0.07) and af ^ o.56 ± 0.31 (0.47 ± 0.41) extracted from ^ D'-*^+p-, and 

B- D(*)V", respectively. Similarly, we use af^ = 0.935 (0.803) and af^ = 0.5 (0.553) 
for the MS form factor case. The CKM matrix elements V^d and Vet are taken to be 0.975 
and 0.039, respectively. 

It is useful to give first an outline of our results. Under factorization, the theoretical input 
for the current-produced D^K^K^, D*^K^K^ modes are all determined, such as from 
EM data and ai from the B^ D^p^ , D*^p^ rates. The calculated rates turn out to be 
in good agreement with the experimental results |l]|. The decay spectra are predictions, 
which can be checked soon. For the D^*^^ K~ case, since we have to fit the unknown 
transition form factors from the corresponding decay rates, it is not as solid as the purely 
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TABLE I: B{B~ D^*^+K-K^) in units of IQ-"^, where the upper and lower Hmits come from 
scanning the Xmin ^ region for the kaon weak form factor . 







MS 


BSW 


Experiment Q 


W- 


D+K-K^ 


-0.21 


-'■•^*-0.20 


1.6 ±0.8 ±0.3 
< 3.1 (90% CL) 


W- 


D*+K-K^ 


9 Q+0.30 
^•°-0.36 


o r)t:+0.32 


2.0 ±1.5 ±0.4 
< 4.7 (90% CL) 



current-produced case. But it is interesting that the resulting decay spectra agree well with 
experimental results We comment on D^*^ K~ K*^ modes before we end this section. 



A. S° D(*^+K-K° 

From Eqs. (Q) and (^, the D^K~K'^ decay amplitude is given by 

A{D^K-K') = ^ V,,V:, ar{pB + Pd) ■ (P2 - Pi) F^i^q') Fl'^iq'), (15) 

where pi, p2 stand for the momenta of K~ and K^, respectively. Factorization implies that 
the amplitude involves only the known weak kaon form factor F^^ and the B ^ D form 
factor F^^ . For 5° D*^K~K^, the decay amplitude is 

A{D*^K-K') = ^V^V:, a,F^^{q') \ ^ ^^^'^ 6,..^ [p, - Pif e^^.p^sP'o* 



+ i 



[niB + mD*)Ai{q ) Ej^* ■ {p2 - Pi) - 



X (pb+Pd*) ■ {P2-Pi] 



(16) 



which involves the B D* transition form factors V, A^ ^- There is no tunable parameters 
for these two modes, hence they provide a good test for the factorization approach. The 
P2 —pi factor implies that the K~K^ pair must be in a P-wave state, which is just what the 
Belle experiment observes 

We show in Table H the calculated rates of these two modes. Both turn out to be at 
the 10~^ level. The upper and lower limits are from scanning the Xmm + ^ region for the 
maximum and minimum of the kaon weak form factor F^^. We find good agreement 
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FIG. 4: The K~K^ mass spectrum for 5° D~^K~K^ (lower) and D*^ K~ K'^ (upper), where 
sohd (dashed) hue stands for using the MS (BSW) hadronic form factors. 

between factorization and experimental results, especially for the D^K~K^ case. This 
provides evidence that the factorization approach works. 

It is instructive to understand the contribution to the D^*^^ K~ K'^ rates from the resonant 
and non-resonant parts of Ff^, where the latter refers to the Xi^2 and yi terms in Eqs. (^ 
and (H). As previously noted, the p contribution dominates the resonant part. We find that 
the resonant part contributes 40% (43%) of the D^*^^ K'^ rate, while the non-resonant 
part contributes 13% (15%). Constructive interference between the two is needed to give 
rates close to experimental results, hence is indeed important. 

The K~K^ mass spectra of the D^*^'^ modes are shown in Fig. |^. Both peak 
close to threshold, which is due to the near-threshold behavior of the F^^ form factor (see 
Fig. ^). There is no other clear structure, other than the B D* form factor effect at larger 
q^. Because of lower reconstruction efficiencies, the spectra has yet to be measured 

experimentally but our predicted spectrum can be checked soon with more data. 
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B. B- D^*^^K-K^ 



From Eqs. (0), (H) and (W), we have the amphtudes 



02 fow^iq )pb ■ {P2-Pi) 



(17) 



A{D*'K-K') = ^V,,V:,{ a. Frig'' 



^^^^ e,,^p e*j!;,p%pf'^. {P2 -piT 



rriB + mo* 



M\Ai{q^){rnB + rnD,)e*jj, ■ {p2 - Pi) 

\ ^D* ■q{PB+ PD*) ■ [P2 - Pi, 

rriB + ml, 
+02/^ niD* I i w_ (g^) e}), ■ {p2 - Pi) 



+Hq') e,,^f, e*j,^.pUP2 + PinP2 - Pif . (18) 




The ai and 02 terms correspond to current-produced and transition parts, respectively. 
There are only two form factors w_(g^) and /i(q'^) in the transition matrix element because 
K~K^ is seen only in 1~ state. For the D^K~K^ case, only W-{q^) contributes. 

Eq. (|l^) involves only one free parameter in w^{t) = Cyj/t^, which can be obtained 
by fitting the central value of the observed rate B{D^K~K^) = (5.5 ± 1.4 ± 0.8) x 10"'' [g. 
We find c!l?^^^^^^ = -35.4 (-33.0) GeV^ and 109.2 (97.4) GeV^ depending on constructive 
or destructive interference between the current-produced and transition amplitudes, respec- 
tively. If we take /i = for now, Eq. (|T8D would give SMS{bsw) ^ ^ 27 (1.45) x 10"^ and 
35.62 (24.97) x 10"'^. The latter result seems too large compared with the experimental 
result of B{D*^K-K^) = (5.2 ± 2.7 ± 1.2) x 10-^ hence disfavor the destructive Ff^-u_ 
interference case unless a fine-tuned h{t) term is used. To obtain the central value of exper- 
imental D*°K-K^ rate, we find cJ;^s(bsw) ^ 3 ^^^^^-^ ^^y^ ^ (-18.5) GeV^. We 

summarize in Table |T| the relevant parameters and the measured rates by experiment. 

Both D^K~K^ and D*^K~K^ rates are used as input and hence are not predictions. 
However, we can still make predictions on the decay spectra. By taking the fitted values 
of cj^^^^^^'', we plot the differential decay rate for the D^K^K^ mode in Fig. ^ and 

compare with the experimental data. The agreement is good. We see that the data itself 

12 



TABLE II: Fitted values of transition form factor parameters c^j^ and Ch, in units of GeV'^, by 
using the central values of D^*^^ K~ rates. 







MS(BSW) 
Cui_ 


MS(BSW) 


5(10-4) g 


B- - 




-35.4 (-33.0) 




5.5 ± 1.4 ±0.8 


B - 


D*^K-K^ 


-35.4 (-33.6) 


11.3 (13.1) or -16.1 (-18.5) 


5.2 ±2.7 ±1.2 



shows a maximum near the K~K^ threshold, which can be naturally explained by our model, 
where threshold enhancement is a genuine result from the form factors in both current- 
produced and transition processes. If threshold enhancement is even more pronounced, as 
data seem to suggest some structure around mx-K° ~ 1-4 GeV, then perhaps the simple 
approximation of Ci^-(t), h{t) oc 1/t^ has to be reexamined. 

The mass spectrum for D*^K~K'^ mode is plotted in Fig. |] with Ch as given in 

Table ||. This indicates the effect of the additional h{t) term in the transition amplitude. 
Like —>■ D^K~K^, one also has threshold enhancement, which again is a genuine result 
of both the kaon weak form factor and the B^ K^K^ transition form factor. We see 
clearly that the W-{t) and h{t) form factors contribute much in the low K^K^ mass range. 



> 

o 





. 5 



Mk-k" (GeV) 

FIG. 5: B^ D^K^K^ spectrum, where solid (dashed) line is for the MS (BSW) model, and the 
data is from Ref. |l|]. 
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1 1.5 2 2.5 3 

Mk-kC' (GeV) 

FIG. 6: The K~ mass spectrum for — > D*^K~K^, where sohd, dot-dashed, dashed 
and dotted hnes are for MS model with c/j = 11.3, —16.1 GeV'^ and BSW model with Ch = 
13.1, -18.5 GeV^, respectively. 

C. B ^ D^*^K-K*° 

Before we end this section, we comment on some features of B ^ D^*^ K*^ modes 
under the factorization picture. Since there is no independent data such as EM form factors 
that one could use, we do not go into the details. 

The Belle experiment already makes some interesting observations [Q: (i) although 
K~K*^ could have = 0~, 1~, 1+ quantum numbers, the data prefers the = case, 
(ii) the rates seem to be dominated by the ai resonance, but the fitted ai — ^ K~K*^ rate 
is 2-5 times larger than the CLEO result [|l^. The data therefore suggest that the am- 
plitude contains both ai —>■ K^K*^ resonance plus a non-resonant part. The situation is 
quite similar to the D^*^^ K'^ case, where the kaon pair mainly comes from a p-pole plus 
a QCD motivated contribution. The p pole alone only gives about 40% of the rate, while 
constructive interference with the QCD part is quite important. If the same picture holds in 
the D^*^+K'K*^ case, the discrepancy with the CLEO ai ^ K K*^ rate may be resolved. 

On the other hand, the K~K*^ pair is produced by an axial current, and no EM data 
could be used as in the K~K^ case. However, further theoretical tools, such as Weinberg 



sum rules ||2^, may be useful to transfer information on vector current form factors to 
axial current form factor. One can in turn extract axial current form factors from the 
;bo ^ D^*)+K-K*'^ data, which may not be obtained by other means. 
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V. DISCUSSION AND CONCLUSION 



In this paper, we use factorization approach to study three-body B D^*^K~K^ decays. 
There are two mechanisms of kaon pair production, namely current-produced and transition. 
The D^K^K^ and D*^K~K^ modes involve only current-produced contributions. Under 
factorization, the kaon pair can only be generated through weak vector current, which can 
be related to EM current through isospin. These modes provide good means to test factor- 
ization, and the result is encouraging. The D^K~K^ and D*^K~K^ rates in factorization 
approach are in good agreement with data. 

The B- D^*'^^K-K^ decays also receive the transition contribution. The form of 
these transition form factors are determined through QCD counting rules, and we fix the 
parameters by using the measured D^K'K^ and D*^K~K^ decay rates. The predicted 
mass spectra of the B^ D^*^^K^ modes agree well with data and exhibit threshold 
enhancement as do the B^ D^*^^ K~ cases. The D'^K'K^ spectrum shows a peak 
around ttix-ko ~ 1-2 GeV, while D^K^K^ spectrum has a peak around rriK-KO ~ 1-5 GeV. 
This is due to the different l/m^_^o behavior of current-produced and transition form 
factors and can be understood from QCD counting rules. To be specific, QCD counting 
rules require the transition form factors to have a faster damping in ttik-ko spectrum than 
the current-produced one. 

Despite the success in describing the mass spectrum of the D^K^K^ mode, our treatment 
of the B~ — i> K~K^ transition form factors may be oversimplified. Assuming the asymptotic 
form required by PQCD may be too strong an assumption, and might have over-enhanced 
the contribution from the near-threshold region. More careful study on other possibilities, 
such as using pole models for transition via resonances, would be helpful in clarifying the 
underlying dynamics of B^ D^*^'^ transitions. 

It is interesting to compare with the treatment of B^ D*~pn [Q, p^{'K^)pn and 



B^ — s> ppK~ modes |jT8[. The success in explaining the rates of both D*~pn and ppK~ 
modes, and the mass spectrum of the latter, encouraged us to apply the factorization ap- 
proach further to the present cases. However, due to the lack of independent information 
on axial form factors, the difficulty to measure neutrons, and the presence of large number 
of operators in the latter case, none of these modes provide a good testing ground for the 
factorization approach. The present D^*^K~K'^ cases turn out to be a good place for such 
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a test. In particular, the D^K~K^ and D*^K~K^ modes are free from any undetermined 
parameters in the factorization approach. The good agreement with data provide support 
for the idea of factorization plus usage of isospin-related form factor data. Further support 
comes from the agreement of D^K^K^ spectra with data. 

Even a hand-waving physical argument may be welcome to explain why a simple factor- 
ization approach works so well in these potentially complicated three-body decays. Before 
we end this paper, we would like to offer one such argument. The QCD counting rules 
constrain the kaon pair in the light pair mass region. With a small invariant mass, the 
two kaons move colinearly and energetically. This is certainly a conducive situation for the 
kaon pair to decouple from the recoil -D^*-* meson, hence "factorize". This heuristic picture 



therefore resembles the B decay to two meson case |21|. It should be noted that PQCD 



may not be the only explanation of factorization p2| , |23| . The success of factorization in 
5 -> D'-*^K~K° decays urges a serious study of the underlying mechanism. 
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